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Abstract 

We study the kinetics of competitive random sequential adsorption (RSA) of particles of binary mixture of 
points and fixed-sized particles within the mean-field approach. The present work is a generalization of the 
random car parking problem in the sense that it considers the case when either a car of fixed size is parked 
with probability q or the parking space is partitioned into two smaller spaces with probability (1 — q) at each 
time event. This allows an interesting interplay between the classical RSA problem at one extreme (g = 1), 
and the kinetics of fragmentation processes at the other extreme (g = 0). We present exact analytical results 
for coverage for a whole range of q values, and physical explanations are given for different aspects of the 
problem. In addition, a comprehensive account of the scaling theory, emphasizing on dimensional analysis, 
is presented, and the exact expression for the scaling function and exponents are obtained. 

PACS number(s): 05.20.Dd,02.50.-r,05.40-y 

The kinetics of adsorption or deposition of particles on a substrate is one of the most common phenomena 
that occurs in many branches of science and technology including physics, chemistry, biophysics and medicine. 
Examples include adsorption of macromolecules and microscopic particles such as polymers, colloid, bacteria, 
protein or latex particles |^ on solid surfaces. Due to its importance, it has been studied extensively through 
all avenues of research comprising experimental, numerical and analytical means. Owing to the complex nature 
of the process, one can hardly make any progress analytically in more than one dimension and therefore most 
of the analytical work remains confined within one dimension only, where it is exactly solvable. The 1 — d 
problem is not at all a pedagogical matter, instead it has been very useful to provide a deep understanding 
of the underlying mechanism and it has offered an insight into more realistic models in higher dimensions. It 
has also been used as a test model to verify the results obtained from numerical simulations and data achieved 
through real life experiments (for extensive review see |4|). 

The simplest model one can think of that can capture the generic feature of the process is the kinetics of random 
sequential adsorption (RSA) of particles of a fixed size on a clean substrate This is popularly known as the 
random car parking problem. This has both an instructional value as well as a long historical significance. The 
basic rules of the processes are (i) at each time event one particle is adsorbed, (ii) the sites for deposition of 
adsorbing particles are chosen randomly and (iii), once the site is chosen, the event is successful if the site is 
empty; otherwise, the particle goes back to the particle reservoir and competes again for further pick up on an 
equal footing with the rest of the particles in the reservoir. Within this model, no two particles or any of their 
parts can occupy the same spatial position in the substrate. This means that particle overlapping is forbidden 
and hence the resulting model describes the kinetics and structure formation of a monolayer. Moreover, once 
deposited, particles are clamped permanently, making them immobile; therefore, the state of the site at each 
time step is changed irreversibly from an empty to a filled state. The irreversible nature of the process causes 
the lack of a detailed balance and drives the system out of equilibrium. Therefore, one can no longer use the well 
developed standard statistical mechanics approach fit for its equilibrium counterpart. A unique yet fascinating 
feature of the RSA processes is the existence of a jamming limit in the arbitrary dimension. This is the state 
when it is impossible to deposit even one single further particle on the substrate, despite the presence of gaps 
(voids), which are in fact not large enough to allow any further deposition. 

To bring the RSA problem closer to real life experiment, there exist many interesting variants of this model. 
Examples include cooperative sequential adsorption (CSA) ^, the accelerated random sequential adsorption 
(ARSA) model ||7||, the balhstic deposition (BD) model, the RSA on disordered substrates the RSA of 
growing objects [ |10| and many attempts to include transport of depositing particles by diffusion Q. Recently, 



there has been an increasing interest in the study of the RSA of mixture of particles of different degree, 
comprising binary mixture of greatly differing diameters jl^, |l^, and continuous mixture of sizes obeying a 
power-law size distribution fl^ . It is important to note that in the former two cases, the ultimate structure in 
the long time limit is described by the jamming coverage, whereas in the latter case, the resulting monolayer 



is uniquely characterised by an exponent called fractal dimension, since the arising pattern is a fractal |13 
The motivation of our present work comes from Bertelt and Privman who studied the binary mixture of 
monomers and k-mers on 1-d lattice as well as Talbot and Schaaf |l2[ , who studied the RSA of binary mixture of 
greately differing sizes. Bertelt and Privman obtained an approximate expression for the coverage for continuum 
model by taking k — > oo. 

In the present work, we concentrate on the binary mixture of fixed size and point particles, instead of pointlike 
particles, in a 1 — d substrate. We give an exact solution for the gap size distribution function and a continuous 
spectrum of jamming coverage depending on the probability at which cars are parked. It is note worthy to 
mention that within the RSA model only one particle can be adsorbed at each time step. This is the basic 
principle of the RSA process and it is indispensable if we want to call it sequential adsorption. In our investigation 
we consider the case when either a car is parked with probability q or the parking space is partitioned with 
probability (1 — q). That is, points and fixed sized particles compete for deposition, and at each time step only 
one of the two events is successful, and hence we call it competitive sequential adsorption. Indeed, we find 
that both the resulting dynamics and the final coverage differ significantly from previous studies. We find that 
the RSA of mixture of point-like and fixed size particles, points and fixed size particles and binary mixture of 
significantly differing size behave in a completely different way. This means that one has to be very careful in 
making any approximation, since the final results are very prone to the exact ratio of the particle sizes. The 
beauty of the present work is that the results are exact and obtained analytically. 

In addition, we present an explicit scaling description with special emphasis on the dimensional analysis. It is 
worth to mention that, despite the long history of the RSA problem, a proper description of the scaling theory 
has remained untouched, although it has been a potential candidate for this. Dimensional analysis provides an 
insight into the problem and of course one cannot get things wrong if one remains faithful to the dimensional 
consistency. 

We first give the general rate equation that describes the sequential deposition of particles, whose sizes and 
positions are determined by the specific choice of the model. Let us denote P{x, t) as the gap size distribution 
function of size x at time t, consequently t) obeys the following integro-differential equation 

= -Pi^^t) dzp{z) dyF{x-y-z,y\z) + 2 J dyP{y,t) dzp{z)F{x,y - x - z\z). (1) 

Here p{z) is called the parking distribution function, that determines the size of the depositing particle and 
F{x — y — z, y\z) is the deposition kernel, that determines the rate a{x, z) = dzp{z) ^ dyF{x — y — z, y\z) at 
which an interval of size x is destroyed by the deposition of a particle of size z, thus creating two new intervals 
of size X — y — z and y. The first term on the right hand side of the Eq. (1) describes the destruction of the 
interval of size x and the second term their creation from the interval y {y > x). Note that if one chooses 
F{x,y\z) = 1 and p{x) = 6{x — a) {a is the size of the depositing particle) then the Eq. (1) simply describes 
the classical car parking problem [|[ It is further interesting to note that if one chooses p{z) — 6{z) then 
the resulting equation describes the random sequential deposition of points, which is equivalent to the standard 
binary fragmentation of particles p5[ ]. In fact, the car parking problem can be interpreted as placing cuts with 
finite thickness and fragmenting particles into two pieces at each time step. The connection between the RSA 
and the binary fragmentation process was first noted by Ziff p^ . 

In order to gain a detailed insight into the process, we first study the RSA of point-sized particles alone, which 
in fact describes the kinetics of the binary fragmentation process. Let us consider the simplest case when 
F{x,y) = 1, which is known as the random scission model 16 ; the resulting equation is 



^I^^^^xP{x,t) + 2 I P{y,t)dy. (2) 



First, note that a{x) — x is the quantity that describes the rate and therefore it must bear the dimension inverse 
of time t. This means that x and t are interlocked and any one of these two variables can be expressed in terms 
of the other one. On the other hand, since x and t are the only two governing parameters, the quantity P{x,t) 



can also be expressed in terms of one parameter alone. Therefore, assuming t to be the independent parameter, 
we can write the following scaling ansatz 



P{x,t)^Pom{x/s{t)) (3) 

where (/)(^) is known as the scaling function. The Eq. (2) will admit scaling if we can choose a time dependent 
scale s{t) for the spatial variable, and Po{t) for the particle size distribution function so that all plots of 
P{x,t)/ Po{t) against x/s(t) for any initial distribution collapse onto one single curve. This can only be true 
if both Po{t) and s{t) show a powerlaw behaviour. This is due to the fact that Po{t) and s{t) must bear the 
dimensions of P{x, t) and x respectively. On the other hand, the dimension of any physical quantity must be 
of power-monomial nature, i.e. we can write 

Po{t) ~ and s{t) ~ t\ (4) 

Simple dimensional analysis of Eq. (1) immediately reveals z = 1 and since it describes the conservation of 
mass, i.e xP{x, t)dx, it is independent of time and yields a = 2. The solution of the Eq. (1) is now well 
known and reads as 

P(a;,t) Rit^e-^* (5) 

which gives (/)(^) = e^"'. Therefore, both the exact solution and the dimensional analysis confirm that bears 
the dimension of P{x, t) and x bears the dimension inverse of time. We shall use this result later. 

We now study the random sequential deposition of a mixture of points and of particles of size a. Considering 
the rate equation for the classical car parking problem and the binary fragmentation model, we can write the 
following equation 



dP{x,t) 
dt 



= -{x- q(T)P{x, t)+2q dyP{y, t) + 2(1 - q) P{y, t)dy for x>a 

J X-\-(7 Jx 



(6) 



and 



dP{x, t) 
di 



2q I dyP{y, t) + {l-q)(- xP{x, t) + 2 [ P{y, t)dy) for 



x<a. (7) 



This is an obvious generalization of the RSA problem, in which at each time event a particle of size <j is deposited 
with probability q or a point is deposited with probability (1 — q). Another way of interpreting the present 
model is that at each time event a car is parked with probability q on the one dimensional parking space or the 
parking space is divided into smaller spaces with probability (1 — q). Knowing the solution of Eq. (2) and the 
solution for the car parking problem alone, we write the following ansatz 

P(a;,t) = A(i)e-(^-«'^)-^W (8) 

where A{t) and B{t) are yet to be determined. The simplest way of finding B{t) is by appreciating that the 
argument of an exponential function must be a dimensionlcss quantity. On the other hand Eq. (6) implies that 
{x — qa) is the rate at which gaps of size x are destroyed, and hence it must bear the dimension inverse of 
time. Therefore, only B{t) ~ t can make the solution physically acceptable. In order to obtain the solutions for 
A{t) and B{t) one can also substitute the ansatz into the Eq. (6). This gives two differential equations, (i) the 
differential for B{t), ^^^^ = 1 and (ii) the differential equation for A(t) 

d\nA{t) „ e""^* 

^ = 2,^ +2(1-,) A. (9) 

Both the equations should be solved subject to the initial condition that we start the process with an empty 
substrate i.e. P{x, 0) = 0. Therefore, the solution for B{t) is trivial, which only confirms the solution we 
obtained through the dimensional analysis. Note that the solution of the equation for A{f^ must be of the form 
A{t) oc t^, since P{x, t) must bear the dimension of t^ . It is then quite straightforward to obtain the solution of 
Eq. (9) which is 

A{t) = t''Fg{at) (10) 
where Fg{at) = e 'Jo " Therefore, the solutions of Eqs. (6) and (7) are 

P{x,t) =t^Fq{at)e-^''-i'''^' for x>a (11) 



and ^ 

P{x, t)= [ dTTFq{aT)e~^''^'"'^^ bqe-"^ + (2 - t)(1 - q)] for x < a. (12) 
Jo L 

The above solutions imply that we can still recover the solution of the binary fragmentation process by taking 
(7 = 0. It is important to mention that the solution of Eq. (6) alone is enough to provide us with all the 
interesting information we need. In fact, we never use the solution for P{x, t) when x < a, which is also true in 
the case of the classical car parking problem. 

To give a scaling analysis and for better clarity we define P{x, t) = P{{x — qa), t) for x > a, then we find that 
the solution of Eq. (6) satisfies the following identity 

P{{x - qa)X, t/X) = \-'^P{{x - q(j),t). (13) 

This relation is the hallmark for the existence of scale invariance and it is equivalent to a data- collapse formalism 
JlSt , as we shall show below. It means that if the deposition rate (x — qa) is increased by a factor A, and the 
observation time is decreased by the same factor, then the resulting structure would look the same except for 
a numerical prefactor. In other words, increasing {x — qa) by a factor of A and decreasing t by the same factor 
means that the numerical value of the gap size distribution function is reduced by a factor of A^ . Since the Eq. 
(13) is true for all positive values of A, we can choose A = i in Eq. (13) then 

Pscale^P{Ll)^<l>[Ul) (14) 

where 

P 

Pscaie = ^; ?? = 0-i and ^ = (x - qa)t (15) 

and (/)(^, 77) is the scaling function. Comparing with the explicit result we get 

cl^{^,v) ^ Fg{r^)e~^ (16) 

and hence (j) is indeed a dimensionless quantity. Therefore, if Pscaie is plotted against the scaled size ^, then 
the entire family of curves collapses into one single curve described by the function ^(^,77). This is possible 
due to the fact that the scaling function depends only on a dimensionless quantity, the numerical value of 
which must be independent of the choice of units. The method of data collapse is a very powerful technique 
for establishing scaling and it is especially useful to analyze and extract exponents from data obtained from 
numerical simulations or from real life experiments. 

The most important quantity of interest in the RSA process is the coverage or jamming limit which is defined 
as 

/>oo 

9{t) = 1 - / xP{x,t)dx. (17) 
Jo 

It tells us the fraction of the substrate covered by the depositing particles. In the case of classical random car 
parking problem only 74.759 percent of the total substrate is covered by the depositing particles. It is now 
interesting to see, how the deposition of points or the division of parking spaces changes the jamming limit. In 
the present case, of course, there is always space for adsorption of points as the points do not have any width 
due to their inherent definition. Hence, they do not contribute to the coverage or the jamming limit. Therefore, 
in this model the jamming limit is the state when there are no more gaps available for the adsorption of particles 
of size a. Once this state is achieved, it is no longer necessary to continue the process as far as the jamming 
limit is concerned. So, it would be interesting to see in what way the coverage changes, if at all. 

To find this out, it is more convenient to deal with the rate equation for the coverage than the coverage itself. 
Combining Eq. (6-7) and Eq. (17) we get 

d9 

— = qa dx{x - a)P{x,t). (18) 
o-t J ^ 

Clearly it states that only fixed size particles contribute to the coverage. Notice that the right hand side of Eq. 
(18) bears the dimension inverse of time, as it should be. Substituting the solution of Eq. (6) into this and 
integrating it, we find the final expression for the jamming limit 



/•oo 

e{oo) = q dsFg(s)e-^^-i'>\ 
Jo 



(19) 



This is a dimensionless quantity and hence, upon transition from one unit of measurement to another within 
a given class, its numerical value must remain unchanged. Therefore, the final coverage or the jamming limit 
6{oo) is independent of the size of the depositing particle cr, which is indeed a non-trivial and interesting result. 
This also reflects the fact that the gap size distribution function P{{x ~ qa), t) satisfies the exact identity (Eq. 
(f3)) and the function Fq{at) is a dimensionless quantity. In fact, all properties of interest are independent of 
the size of the depositing particle, and the rate at which they are adsorbed is due to the existence of the scaling 
property. Clearly, if a is increased by a factor A, then the time to reach the jamming configuration is reduced 
by a factor A. However, this is true only if the substrate is sufficiently large in comparison to the depositing 
particle. One can now obtain a whole range of values for the jamming limit by just tuning the q value. Fig. 1. 
shows that the jamming limit starts from zero at g = 0, since there the problem reduces itself to the random 
sequential adsorption of points. Note the inherent definition of point that does not occupy space. Therefore, 
one might apprehend that it does not play any role in achieving the jamming limit if it is adsorbed sequentially 
with a fixed size particle. However, note that once a point is deposited, it excludes fixed-sized particles from 
landing within around cr/2, which is the distance between the point where it is chosen to land and the centre 
of the fixed sized particle. This exclusion will certainly have effect on the resulting dynamics of the process. 
In other words, a point divides the substrate, and hence each successful deposition of points creates two new 
ends belonging to two different gaps. No particle can be adsorbed on a space which belongs to two different 
gaps i.e. overlapping of points is forbidden. We find that for q > 0, the jamming limit changes nontrivially. 
Fig. (1) exhibits a 3 dimensional view of the jamming limit against time and the probability q with which 
fixed-size particles are adsorbed. Fig. (2) is the cross-section of the coverage-time plane at g = 1 and q = 1/2 
of fig. (1). It is generally believed that if a mixture contains a small number of particles of different sizes, 
then the kinetics and jamming limit are primarily determined by the smallest size which in the present case is 
a point particle. However, in our study we find that the jamming limit is primarily determined by the larger 
particle; yet its dynamics and other aspects are influenced by the size of the smaller particle. For q > 0, the 
jamming limit is always smaller than its corresponding classical counterpart. Fig. (2) is plotted to demonstrate 
the slower approach to the jamming limit as q value decreases. This is due to the fact that the adsorption 
of points leads to the crowding of more partition in the parking space than to the adsorption of fixed sized 
particles. Fig. (3) in fact is the cross-section of coverage-g plane of Fig. (1) at the state when the system has 
achieved the jamming limit. It shows that the jamming limmit increases monotonously as fixed-size particles 
wins more in the competition with the point particles for adsorption. However, it does not increase linearly 
with q, instead one can subdivide the whole curve into a number of regimes of width 0.2 in q values, so that 
in each regime they can be well approximated to a straight line. However, notice that, between q = 0.95 and 
g = 1, the coverage increases much stronger than in any previous regime. The appearance of a non- linear rise 
becomes more pronounced if it is plotted in the larger scale. 

In conclusion, the present model, though simple, can yet capture some generic features of the RSA of mixture 
of particles such as coverage and the scaling behaviour. We find that the coverage of mixture of points with 
fixed-sized particles stays always lower than its classical counterpart. It would be interesting to see how the 
results change when both the contents of the binary mixture are of a finite size but differ in length which we 
intend to present in a forthcoming paper [ p^ . In addition, we give an extensive scaling description of the RSA 
problem which has not been addressed so far. One potential application of the present work could be the simple 
cyclization reaction process where adjacent pendant groups link randomly along the polymer, or the polymer 
itself undergoes a possible degradation. This is in fact a natural generalization of the original work by Folry 
1^ . We believe the present work will shed a new insight into the underlying mechanism of the problem and will 
help guiding numerical and experimental works dealing with more realistic situations addressing the adsorption 
of mixture. 
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grateful to the Alexander von Humboldt Foundation for awarding the fellowship. 
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Figure 1: Jamming limit as a function of time t and probability q (a 3D view). 
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Figure 2: Jamming limit 6{qo) against time t. The lower curve is for q ~ ^ and the upper curve for q — 1 (the 
classical car parking). 




Figure 3: Jamming limit 6'(oo) as a fmiction of q (this is in fact [100] face of the Fig. (1)). 



